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Measurement of Rb 5P 3 / 2 scalar and tensor polarizabilities in a 1064 nm light field 
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We employ doubly-resonant two-photon excitation into the 74S Rydberg state to spectroscopically 
measure the dynamic scalar polarizability, Qo, and tensor polarizability, ay, of rubidium 5 P 3 / 2 . To 
reach the necessary high intensities, we employ a cavity-generated 1064 nm optical-lattice light 
field, allowing us to obtain intensities near 2 x 10 11 W/m 2 . In the evaluation of the data we 
use a self-referencing method that renders the polarizability measurement largely free from the 
intensity calibration of the laser light field. We obtain experimental values ao = —1149 (±2.5%) 
and «2 = 563 (±4.2%), in atomic units. Methods and results are supported by simulations. 

PACS numbers: 32.80.Rm, 37.10.Jk, 32.10.Dk, 37.30.±i 


Polarizabilities of atomic energy levels govern the re¬ 
sponse of an atom to an external electric field and are a 
fundamental property to be considered in atom-trapping 
and precision-measurement experiments. For example, 
the polarizabilities are required for determining magic 
wavelengths of state-insensitive trapping in optical lat¬ 
tices [1, 2]. Theoretical calculations of dynamic polariz¬ 
abilities are complicated, and yet available experimental 
measurements might carry large uncertainties due to the 
difficulty in calibrating the field strength experienced by 
the atoms. Here, we report a measurement of rubidium 
scalar and tensor polarizabilities conducted in a strong 
1064 nm light held, where the magnetic sublevels of Rb 
5 P 3/2 are resolved and the data analysis is largely free 
from the calibration of laser intensity. Despite the fact 
that 1064 nm optical traps for rubidium atoms are widely 
used in cold-atom trapping, to our best knowledge there 
is no other such experimental measurement. Our work 
is not only applicable to experiments utilizing Rb 5P 3 / 2 
levels in 1064 nm laser traps, but also serves as an exper¬ 
imental test for validating and improving existing theo¬ 
retical models for the polarizability. 

In the presence of an optical held, an atom is polarized 
and its energy levels are shifted. The applicable atom¬ 
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FIG. 1. (color online). Experimental details, (a) A near- 
concentric cavity, which has a focus at its cavity center, is 
used to generate a GHz-deep 1064 nm optical-lattice laser 
trap, (b) Timing sequence showing lattice intensity, Rydberg- 
atom excitation and field ionization pulses vs time (t = 30/rs, 
repetition rate = 100 Hz), (c) Key components of the PDH 
cavity stabilization. 


held interaction Hamiltonian, He, in a linearly-polarized 
electric held with amplitude E 0 is 


B 0 I T , , T if / \ , / \ — J(J + 1) 

X \J, r n J )(J,mj\\a 0 {u}) + a 2 {u}) — j( 2 J -1) 


( 1 ) 


where ao(u>) and < 22 (to) are frequency-dependent a.c. 
scalar and tensor polarizabilities, and J and mj are quan¬ 
tum numbers of the total electronic angular momentum. 
The dynamic polarizability of the rubidium ground-state 
5S!/2 depends only on ao (a 2 is zero because J = ^). For 
5 P 3 / 2 , both ao and a 2 contribute to the polarizability. 
The full Hamiltonian includes the hyperhne structure, 
H = H H fs + He, with 


Bhfs= 4Lhfs I • J ± 

3(i-J) 2 ±|l-J-/(/±l)J(J±l) 
HFS 2/J(2/ — 1)(2J — 1) 


( 2 ) 


where I is the nuclear spin, Ahfs is the magnetic-dipole 
and .Bhfs the electric-quadrupole hyperhne constant. 
Hhfs and -Bhfs are well known for 87 Rb and 85 Rb [3, 4]. 
The octupole contribution is omitted because it is too 
small to be observed here. 

We measure ao and a 2 using Rydberg two-photon exci¬ 
tation spectroscopy. The data analysis is based on linear 
hts of spectral data sets in a modihed a.c. Stark map, in 
which the frequencies of the two excitation lasers are plot¬ 
ted against each other, with the unknown polarizabilities 
as htting parameters. Our method has the strength that 
it does not require a precise calibration of the 1064 nm 
laser intensity at the atom trapping site. Also, light shifts 
are in the range of several GHz, which is important for 
a precise measurement of ao and a 2 of Rb 5P 3 / 2 (level 
width 6 MHz). 

We utilize an in-vacuum, moderate-hnesse (« 600), 
near-concentric optical cavity [see Fig. 1 (a)] at 1064 nm 
to generate deep optical-lattice potentials in a linearly- 
polarized light Held (for details see Ref. [5]). We load the 
lattice at low intensity directly from a 87 Rb magneto- 
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optical trap (MOT), and then increase the lattice power 
after the MOT light is pulsed off [timing see Fig. 1 (b)]. 
The smooth sinusoidal lattice-intensity ramp compresses 
the atom distribution in the lattice wells, where light 
shifts reach several GHz. Atoms are excited into Ryd¬ 
berg states while the lattice is held at a high, constant 
intensity. Rydberg atoms are detected by field ioniza¬ 
tion [6] and ion counting. After detection, the lattice is 
switched back to the loading intensity for the next ex¬ 
perimental cycle. During the loading phase, the 5 Si/2 
ground-state atoms have an estimated trap frequency of 
about 700 kHz along Z and about 7 kHz along X and Y. 
Hence, the 30 /is lattice ramp leads to adiabatic compres¬ 
sion in the Z and mixed adiabatic/diabatic compression 
in the transverse directions. 

The in-vacuum near-concentric lattice resonator is sta¬ 
bilized to the 1064 nm trap laser (short-term bandwidth 
100 kHz) by the Pound-Drever-Hall (PDH) scheme [7] 
[see Fig. 1 (c)]. The electro-optic modulator (EOM) gen¬ 
erates PDH frequency side bands. The reflection from the 
cavity is sampled by a beam sampler for synthesizing the 
PDH error signal. The PDH feedback circuit has two out¬ 
puts: a high-voltage slow feedback is applied to a piezo, 
which translates one of the cavity mirrors, and a fast 
feedback frequency-modulates the acousto-optic modula¬ 
tor (AOM), which compensates the rapid frequency fluc¬ 
tuations of the trap laser. The lattice intensity ramp is 
generated by amplitude-modulating the AOM. To ensure 
acceptable performance of the lattice cavity lock dur¬ 
ing lattice compression, the amplitude of the PDH error 



FIG. 2. (color online). Experimental spectra overlaid over 
a simulated spectral-density plot. The grey drop lines mark 
the (Abiue,A re d) coordinates of the peaks of the evaluated 
spectral lines. The bold, short black lines underneath the 
evaluated spectral lines indicate the uncertainties of the 
peak coordinates. The dashed grey lines show linear fits 
through the marked positions (Abiue, A re d). For the simulated 
density plot, the peak lattice intensity is 1.8x 10 11 W/cm 2 , 
Qo(5P 3/2 ) = -1149.3, a 2 (5P 3 /2) = 563.3, and a 0 (5S 1/2 ) = 
687.3 (polarizabilities in atomic units). 


signal is normalized in real time by a trap-laser power 
reference. The normalization circuit is composed of op- 
amp logarithmic and exponential amplifiers. Our method 
to stabilize the optical cavity at high intra-cavity power 
differs from the strong-weak beam method [8] and from 
methods that use a separate reference cavity [9]. 

The polarizabilities are derived from two-photon step¬ 
wise (resonant) excitation signals from the light-shifted 
5 Si/ 2 ground state through a light-shifted 5P 3 / 2 sublevel 
into the 74 Si/2 Rydberg state. For each of the selected 
lattice intensities, the two-photon excitation spectra are 
taken for a set of fixed lower-transition detunings, A rec j, 
by recording Rydberg counts as a function of the upper- 
transition detuning, Abiue- We reference Abiue to a nar¬ 
row spectral line that corresponds to off-resonant two- 
photon excitation of Rydberg atoms outside the lattice. 
This line serves as a well-defined reference point. For 
each lattice intensity, the spectra are arranged in water¬ 
fall plots in which the //-displacement is given by A rec j, 
while Abiue is plotted along the x-axis. In Fig. 2 we 
show such a modified a.c. Stark map measured for a 
transmitted lattice power of 20 mW, which corresponds 
to a peak intracavity intensity at the lattice sites of about 
1.8x 10 11 W/m 2 . The uncertainty of the transmitted lat¬ 
tice power is about 8% and that of the peak intracavity 
intensity is even larger. It is an essential advantage of our 
method that these uncertainties do not affect the values 
of «o and 0:2 obtained from the data. This is because ao 
and «2 are derived from slopes in the modified a.c. Stark 
map, not from absolute level positions. In the following, 
we explain this self-referencing method in detail. 

In Fig. 3 we show the light shifts of the 5 Si/ 2 —>5P 3 / 2 
transitions versus intensity, Iq = ^ceo-Eg, calculated by 
diagonalization of the full Hamiltonian (Eqs.l and 2) us¬ 
ing approximate values of ao and a.^- The ao and a 2 
are the same for both rubidium isotopes, because they 
only depend on electron-field coupling. Since the 5Siy 2 
state has no a.c.-split sublevels, the splitting is only due 
to the 5P 3 /2 state, which has 16 sublevels for s 'Rb and 
24 for 85 Rb. The plots exhibit three intensity regimes. 
In the weak-field regime, Hjjfs dominates, and |E, tuf) 
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FIG. 3. A.c. Stark shifts of the transitions 5 Si/ 2 —> 
5P 3 / 2 of 87 Rb ( 85 Rb) relative to the transition 5 Si/ 2 , F = 
2(3) —> 5P 3 / 2 , F' = 3(4), calculated for cto(5Si/ 2 ) = 700, 
«o(5P 3 / 2 ) = —1100, and a 2 (5P 3 / 2 ) = 550 atomic units. 



















3 


(F = J + I) is the “good” basis. Levels with the same F 
and m/? [-values are degenerate, and their energy shifts 
are linear in intensity. At intermediate intensity neither 
.Hhfs or He dominates, and the energy levels are gener¬ 
ally nonlinear. The level crossing at 7 x 10 10 W/m 2 (for 
8 ' Rb) depends mostly on a 2 and the zero-field hyperfine 
splittings. For the experimental determination of ao and 
a 2 we utilize the strong-field regime, which is analogous 
with the “Paschen-Back” regime of the Zeeman effect. 
With our near-concentric cavity setup, we comfortably 
reach the strong-field regime. 

In the strong-field regime He dominates, and 
| IJmimj) becomes the “good” basis. The energy lev¬ 
els become linear functions of intensity again, and they 
separate into two groups of fixed \mj\. Since a 2 > 0, the 
group with higher energy includes levels with \mj\ = 
while the group with lower energy has |toj| = |. Within 
each subgroup, the energies are split by the residual hy¬ 
perfine perturbation. For example, for 8 ' Rb (7 = | 
and J = |), in the subspace \rnj\ = | the off-diagonal 
terms of I ■ J and (I ■ J) 2 vanish in the |to/TOj) basis, 
and the residual hyperfine-induced shifts are given by 
the diagonal terms (which follow from mjmj = or 
±|). Thus, the lower subgroups in Fig. 3 appear nearly 
equally spaced, due to the leading magnetic-dipole hy¬ 
perfine term Ahfs I • J- The deviation from an equal 
spacing is due to the electric-quadrupole hyperfine term, 
whose strength is several percent of that of the magnetic- 
dipole term. In the subgroups with |mj| = \ the off- 
diagonal elements of I • J are generally non-zero, and the 
residual hyperfine shifts are not approximately equidis¬ 
tant. In our experiment we choose 87 Rb, because the 
number of rn; sublevels is less and the residual hyperfine 
splittings in the Paschen-Back regime are larger than for 
85 Rb. However, the method also applies to 85 Rb. 

The 5S!/2 and 5P 3 / 2 levels exhibit a local response to 
the 1064 nm lattice light field, i.e. their light shifts are 
given by — (a/4)\E 0 (R)\ 2 : where E 0 (R) is the field ampli¬ 
tude at the center-of-mass location J? of the atom. Due 
to their size, Rydberg atoms have a non-local response 



FIG. 4. (color online). Relevant levels of the two-photon 
double-resonant Rydberg excitation, (a) Outside lattice. (b,c) 
Light-shifted levels (solid lines) in transverse (b) and in¬ 
direction (c) in the 1064 nm lattice. The “magic” 74S state 
has a constant light shift along Z. Dotted lines indicate the 
unshifted levels. The intervals between levels are not to scale. 
The arrows illustrate several instances of doubly-resonant ex¬ 
citation on one of the lines in Fig. 2 (see <0> in Fig. 2). 


to the field. The Rydberg-atom light shift is 

Kd(R ) = ~\a e J |£ 0 (R + r)| 2 |V>(r)| 2 d 3 r (3) 

where a e = —545 at. un. follows from the free-electron 
ponderomotive energy [10]. Thus, the Rydberg-atom 
light shift is an average of the free-electron response, 
with the Rydberg electron’s probability density |t/>(r)| 2 
as a weighting factor. While generally t/;(r) also depends 
on R, for the non-degenerate nS Rydberg levels it is 
t/>(r) = (r|nS). The averaging in z-direction is important, 
because the size of the Rydberg atom is on the order of 
the lattice period. The averaging in p direction is not 
important, because the cavity-mode waist Wq is much 
greater than the size of the Rydberg atom. For a few 
“magic” states, such as Rb 74S;iy 2 in a one-dimensional 
1064 nm lattice, I4d only depends on p and not on Z. 
For the magic states, 

Kd (P,Z) = Pad (p) = Y^Lx ex P(^y-) ( 4 ) 

where 7£ max is the peak field amplitude in the entire lat¬ 
tice. In Fig. 4 we show the position dependence of the 
relevant light shifts as a function of atomic center-of-mass 
coordinates X and Z. The peaks in Fig. 2 correspond to 
doubly-resonant excitation of 74 Si/ 2 through one of the 
multiple 5P 3 / 2 sublevels. 

The density plot in the background in Fig. 2 shows 
the result of a simulation. The simulation accounts for 
the center-of-mass thermal distribution of the atoms in 
the optical lattice (temperature T), which causes most of 
the spectral line broadening. Doppler shifts are negligible 
and are ignored. The simulated count rate is an integral 
over the lattice volume that includes the Boltzmann fac¬ 
tor, exp[—a 5S (i?max - E o (P> Z))/(4fc B T)], and the lower- 
transition saturation parameter, 0.5s/(l + s+4( 6 ^ Hz ) 2 )- 
There, s = I/ 7 sat with saturation intensity 7 sat and 
position-dependent intensity 7. Also, A is the detun¬ 
ing of the 780 nm laser from the light-shifted, position- 
dependent 5 Si/ 2 —t5P 3 / 2 transition frequency. The sim¬ 
ulated count rate is summed over all intermediate 5P 3 / 2 
states. Figure 5 shows a compilation of three simulations 
for 8 ' Rb for different peak lattice intensities and the same 
polarizabilities. The slopes in the modified a.c. Stark 
map are insensitive to the peak lattice intensity. The 
peak lattice intensity affects the signal strength and the 
Abiue cutoff, where the doubly-resonant excitation con¬ 
dition is met at the field maxima (which are at locations 
p = 0 and Z = k x 532 nm with integer k) for all 5P 3 / 2 
sublevels. The peak lattice intensity in the experiment 
can be estimated from the Abiue cutoff; this intensity is, 
however, not required to extract the polarizability values. 

As can be seen from experimental and simulated spec¬ 
tra, the lattice-shifted spectral lines are asymmetric and 
triangular. The line shape results from the variation of 
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FIG. 5. (color online). Simulated spectra. We overlay three 
calculations for the indicated lattice intensities. 


the Rydberg light shift along surfaces of fixed intensity; 
the variation is due to the non-local response following 
Eqs. 3 and 4. The line shape is discussed in [14]. We 
use the positions (Abiue, A red ) of the peaks of the spec¬ 
tral lines as primary markers, as shown in Fig. 2. These 
peaks correspond to combinations of intensity I, Abi U e 
and A red for which the double-resonance condition for 
excitating 74 Si/ 2 is met at locations Z = k x 532 nrn 
with integer k [14]. 

In the modified a.c. Stark map, the peak positions 
(Abiue, A red ) satisfy 


y = A red = A 5 S5P = -(ass - & 5 p)Eq(p, Z = 0 ) 

x = Abiue = A74S + A5S5P — A 5 p 

= A 5 S5P + T( a 5P — 0 !'t 4 s)Eq(p, Z = 0 ) = y — ^ - — 

4 «5S — «5P 


where Eo(p, Z = 0) is the maximal lattice field at a dis¬ 
tance p from the axis. The slopes of the levels are inde¬ 
pendent of Eq and are 


d y _ o?5S — Q5P 
da; a 5 s — «74S 


(5) 


In the high-field (Paschen-Back) regime, the a^p of the 
highest energy levels in each subgroup of rrij \ = | and 
\ m j\ = f are, terms ao and « 2 , 

a 5 p(\rrij\ =-) = a 0 - a 2 (6) 

3 

a 5P {\mj\ = -)= a 0 + a 2 (7) 

The differential and average slopes are 


dy 

dx 




1 ( dy 

2 yda; 


I m o\=h 



2 a 2 

CT5S — C*74S 


ass ~ Q0 

QI5S — a 74S 


( 8 ) 

(9) 


The slopes dy/dx are determined by linear fitting of the 
sets of peaks (Abiue, A red ) associated with the respective 


lines in the modified a.c. Stark map. For rrij = | we 
force the y-intercept to zero because this level has a fixed 
slope through all field regimes and passes through the ori¬ 
gin. For rri 3 = | we fit both the slope and the intercept. 
This is because in the weak field regime the rrij \ = \ level 
connects to |F = 3, mp = 0), which has asp = ao — |a 2 . 
Hence, the high-field fit of the rrij = | level yields 
a slight negative y-intercept, which is on the order of 
—30 MHz for a peak lattice intensity of 1.8 x 10 11 W/m 2 . 
We fit the slopes for measurements at several peak lat¬ 
tice intensities and excitation-pulse durations [14]. The 
numbers of slope samples are 5 and 4 for rrij = | and 
\rrij | = respectively. The weighted average slopes are 
summarized in Table I. The ass and a 74 s are, for the 
present purpose, precisely known (ass = 687.3(5) [11] 
and a 74 s = —272.5(5) [10], in atomic units). Eqs. 8 and 
9 then yield ao = —1149 and a 2 = 563 at. un. 

The uncertainties listed in Table I only reflect the 
linear-fitting uncertainty and the reading uncertainties 
associated with the determination of the peak markers 
(Abiue, A red ). The calibration uncertainties for Abiue and 
A red are 1.6% and 0.2%, respectively. Adding all rela¬ 
tive uncertainties in quadrature, the final uncertainties 
are 2.5% for ao and 4.2% for a 2 [14]. Our result is in 
good agreement with theoretical values ao = —1111.62 
and a 2 = 557.31 obtained in Ref. [11]. Also, we have 
adopted the theoretical calculation of ass = 687.3(5) 
in [11], instead of an earlier calculation from [12] and 
an experimental value ass = 769(61) [13]. 

Our method could be adapted to other atoms in deep 
optical lattices and at different wavelengths, such as for 
Cs in a deep 1064 nm lattice or Rb in a high-power 002- 
laser lattice. Instead of using a magic Rydberg state, 
one may utilize low-lying Rydberg states that are suf¬ 
ficiently small that they exhibit an approximately local 
response to the field (given by the free-electron polariz¬ 
ability). A different approach to polarizability measure¬ 
ment in 1064 nm laser fields can be found in [15]. 

The experimental values for the scalar and tensor 
polarizabilities of 5P 3 / 2 are immediately useful in ex¬ 
periments that require on-resonant transitions through 
5 P 3 / 2 , such as two-plroton preparation of lattice-mixed 
hydrogenic Rydberg states in deep 1064 nm lattices [16], 

TABLE I. Experimental 5 P 3/2 scalar and tensor polarizabili¬ 
ties in 1064 nm light field. The uncertainties are reading and 
fitting uncertainties. Polarizabilities are in atomic units. 


parameter 


value 

uncertainty 

Experimental slope \rrij\ 

_ 3 

2 

1.33 

4.19xl0" 3 

Experimental slope rrij 

_ 1 

2 

2.50 

4.47xl0- 2 

differential slope 


1.17 

4.49xl0- 2 

average slope 


1.91 

2.24X10- 2 

ao (experiment) 


-1149 

22 

a .2 (experiment) 


563 

22 
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Rydberg-EIT in 1064 nm optical traps, and evaporative 
cooling utilizing 1064 nm light fields. 
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SI. SPECTRAL LINE SHAPE 

The triangular shape of the lattice spectral lines results 
from the non-locality of the Rydberg atom’s response to 
the Held, which leads to a variation of the Rydberg-level 
light shift along surfaces of constant intensity within the 
optical lattice. This is explained in the following. 

Along the lattice axis (Z- axis) the lattice laser inten¬ 
sity sinusoidally varies with a period of 4 =532 nm. In 
our coordinate system, the lattice-intensity maxima are 
located at Z = k x 2 , with integer k. In the transverse 
direction, i.e., the X- and Y-axes, the lattice intensity 

follows a Gaussian profile lip, 0) = / m axexp(^£-) with 

_ W 0 

a beam waist Wq of 23 pm and p = \JX 2 + Y 2 . On 
any plane that contains the lattice axis, the lattice in¬ 
tensity distribution can be described by equal-intensity 
surfaces that resemble the shape of ellipses, as shown in 
Fig. SI (a). At the locations (p = 0, Z = k |), with inte¬ 
ger k, the lattice laser intensity is maximal, / = J max . 
For decreasing intensity, the equal-intensity “ellipses” 
become larger and increasingly eccentric, as shown in 
Fig. SI (a). 

The atoms on the same equal-intensity surface have the 
same lower-transition detunings and Boltzmann factors; 
however, they experience different lattice-induced light 
shift for the 74S magic Rydberg state. The fight shift of 
the 74S magic state is constant on cylindrical surfaces of 
constant p, i.e. the 74S light shift is given by the lattice 
intensity I(p, 0) = I max exp( = %-). As a result, as shown 
in Fig. Sl(a), on any given constant-intensity “ellipse” 
with I < 7 max there are two extreme cases: point A, at 
the major axis, and point B, at the minor axis. Atoms 
at point A exhibit minimal light shift of the 74S state, 
and the signal strength due to atoms near A is maximal 
because the volume associated with the corresponding 
region is the largest. Atoms at point B experience max¬ 
imal light shift of the 74S state, but the signal strength 
is small. 

Recalling that the 5S and 5 P 3/2 shifts are fixed on 
equal-intensity surfaces, the upper-transition detunings 
on equal-intensity surfaces vary over a frequency range 
L given by the variation of the 74S fight shift, as shown 
in Fig. Sl(b). Figure Sl(b) also shows that the range 
L increases with decreasing value of the lower-transition 
detuning, A 5 g + A 5 p, which is seen as follows. The 
lower-transition light shift, A 5 S + A 5 P, has to approx¬ 
imately equal the 780-nm laser detuning A re( j for one 


of the light-shifted intermediate 5P 3 /2 levels, so that 
level becomes resonantly populated. These resonances 
occur at certain intensities 7 res . The 74S light shift along 
the equal-intensity “ellipse” for 7 res then varies between 

- i a74S jio Imax at point B and ~2 a74S j|7 res at point 

A. Hence, the frequency range L associated with that 
instance of a doubly-resonant Rydberg excitation is 

1 2 

7* = 7^745 [4 max Ires] (SI) 

4 ceo 

This finding accords with the observation that for smaller 
lower-transition detunings A re( j the lattice spectral lines 
become wider. In addition, since the volumes associated 
with points along the equal-intensity surfaces increase in 
proportion with p, the signal increases when moving from 
points B to points A. As a result, with decreasing A red 
the lattice spectral lines become more asymmetric and 
triangular-shaped, as illustrated in Fig. Sl(c). 

Most importantly, the spectral-line markers 

(A re d,Abiue) used in the paper correspond to loca¬ 
tions of the type A, where the lattice spectral lines 
peak. Hence, all relevant light shifts (5S, 5P 3 / 2 , and 
74S) that correspond with the markers (A re d, Abi ue ) are 
proportional to the same intensity I(p, 0). As a result, 
in our self-referencing method the intensity I[p, 0) drops 
out, allowing us to relate the 5P 3 / 2 polarizabilities to 
slopes in the modified a.c. Stark map. This proves one 
of the strongest features of our method, namely that the 
uncertainty of the calibration of lattice intensity against 
optical-lattice power injected into the cavity does not 
affect the uncertainty of the measured polarizability 
values. 


S2. STATISTICAL ERROR 

We conduct the measurement at different lattice inten¬ 
sities and R.ydberg-excitation durations (2 ps or 10 ps), 
as shown in Fig. 2 and Fig. S2. With 10 ps Rydberg- 
excitation pulses, the lattice spectral lines tend to be 
saturated, but this does not shift the locations of the 
peak centers. Each data point (Abi ue , A re d) in the 
modified a.c. Stark map has an individual uncertainty 
(d(Abiue), d(A re d)), which is due to the uncertainty as¬ 
sociated with both excitation-laser frequency stabiliza¬ 
tion methods, the laser line widths, and the uncertain¬ 
ties in locating the peaks of the individual lines. We use 
OriginLab linear regression to determine the slopes b of 
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FIG. SI. (color online), (a) Constant-intensity “ellipses” for intensities 7 max > 1 1 > 7 2 . The eccentricity of the ellipses increases 
with decreasing intensity, (b) Corresponding light shifts of the relevant levels. The shift of the Rydberg state varies over a 
range L that increases with decreasing intensity, (c) Typical triangular line shape. The left width 8 is small and given by the 
level width of the intermediate state (6 MHz) and the laser line widths (~ 1 MHz). The right width is largely given by L. 


the highest energy levels in the \rrij\ = \ groups and the 
|nij | = | groups. The fitting returns a fitting uncertainty 


<5&i. 


The effects of the reading uncertainties 
(d(Abiue),<5(A roc j)) is calculated as follows. Assum¬ 
ing there are M data points for one energy level, 
with x = Abiue and y = A re d- Each data point 
(x = Abiue, y = A re d) carries a Gaussian point spread 
function, with cr x = d(Abiue) and a y = <5(A re d)- We 
generate 10 6 sets of data points obtained with these 
point spread functions, and using random numbers. 
For each set, the slope b of the best-fit linear function, 
y = a + bx , is 
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The 10 6 values of b yield a distribution that has a stan¬ 
dard deviation <5£>2, which constitutes the uncertainty in 
b due to the reading uncertainties (<5(Abi U e), <5(A re d)). 

The net statistical uncertainties 5b of the slopes are 
the quadrature sums of 6bi and d& 2 • In Table SI we 
show V and 5b l for the measurements performed, where 
i is a counter for the measurement sets. We have five 
sets of scans, and obtain TV = 4 slopes for the highest 
level of the group m, = and TV = 5 slopes for the 
group rn-j | = |. The averaged experimental values of 
the slopes and their uncertainties, b and 5b, are then 


calculated using weighted averages, 


rELW 

Eti m ~ 2 


5~b 2 = 


Ef=i m~ 


(57) 

(58) 


In the last step of the uncertainty analysis for the 
slopes, we add the relative uncertainties 5b/b and the 
scaling uncertainties, explained in the next section, in 
quadrature. This yields the net total uncertainties 5b net 
for both slopes. The final uncertainties for the polariz¬ 
abilities are then calculated using standard uncertainty 
propagation in Eqs. 8 and 9. 


S3. SCALING ERROR 

The lower (780 nm) and upper (480 nm) transi¬ 
tion lasers of the Rydberg two-photon excitation are 
frequency-stabilized with the transmission spectra of con- 
focal optical cavities. The lower transition utilizes a con- 
focal cavity at 780 nm, whose free-spectral range (FSR) is 
374.68 MHz with uncertainty of 0.50 MHz. This number 
is calibrated using the frequency reference of the well- 
known 85 Rb transitions: 5 Si/ 2 -A 5P 3 / 2 F = 2 —> F' = 3 
and F = 3 —>• F' = 4. The scaling error of A re( j is less 
than 0.2%. 

The upper transition 480 nm laser is frequency-doubled 
from a 960 nm laser, which is stabilized to a pressure- 
tuned confocal cavity at 960 nm. The FSR of the 960 nm 
confocal cavity is 489.37 MHz with uncertainty of 0.56 
MHz. The calibration is referenced to a RF function 
generator (HP 8656B). The scans of the 480 nm laser are 
done by pressure-tuning the 960 nm FPI, which is driven 
by a step motor. In each of the scans, the raw data are 
recorded in terms of step numbers. We convert the step 
number into a MHz offset at 480 nm; the conversion is 
also referenced to a RF function generator. We found 
that the conversion factor from step number to MHz in 
480 nm is not perfectly linear due to the nonlinearity of 
the pressure response to the step motor’s position. Dur- 
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FIG. S2. (color online). Waterfall plots (modified a.c. Stark maps) of experimental spectra. The grey drop lines mark the 
(Abiue, A re d) coordinates of the peaks of the evaluated spectral lines. The bold, short black crosses underneath the evaluated 
spectral lines indicate the uncertainty on the peak positions in both Abiue and A re d- The dashed grey lines show linear fits 
through the marked positions (Abi ue , A re d). The transmitted lattice powers and the Rydberg excitation pulse durations are 
(a) 32 mW and 10 /zs, (b) 36 mW and 10 /zs, (c) 36 mW and 2 /zs, (d) 40 mW and 10 /zs. 


ing a single scan the nonlinearity is less than about 1.5%. The resultant total scaling error of Abiue is estimated to 

be less than 1.6%. 
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TABLE SI. Slopes of the linear fits in the modified a.c. Stark maps from Fig. 2 and Fig. S2 and their uncertainties. The 


sample number is 

5 for the slope of 

\rrij \ = |, and 4 for the slope of | rrij 

1 = 1- 



figure 

power (mW) 

pulse (fxs) 

slope |toj| = | 

error 

slope \m,j\ = \ 

error 

Fig.2 

20 

2 

2.45 

6.95 x 10 -2 

1.32 

1.05 x 10" 2 

Fig.S2 (a) 

32 

10 

2.58 

9.46 x 10" 2 

1.33 

1.13 x 10" 2 

Fig.S2 (b) 

36 

10 

2.53 

8.15 x 10" 2 

1.33 

1.01 x 10" 2 

Fig.S2 (c) 

36 

2 

- 

- 

1.32 

8.14 x 10" 3 

Fig.S2 (d) 

40 

10 

2.47 

1.78 x 10" 1 

1.33 

8.17 x 10" 3 



